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SUMMARY 

A general derivation is presented for the radiation from a uni- 
formly radiating sphere incident on a small flat plate, spinning about 
an axis not coincident with its normal. The results are presented as a 
function of the separation of the bodies, the orientation of the plate, and 
the orientation of the spin axis. In addition, a series of curves is given 
which represents the power input from earth radiation to one side of a 
spinning flat plate (averaged over a spin period) for a range of altitude 
from 200 to 32000 km. These curves are based upon the assumption 
that the earth is a uniform diffuse emitter radiating as a blackbody at a 
temperature of 250°K. 
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POWER INPUT TO A SMALL FLAT PLATE 
FROM A DIFFUSELY RADIATING SPHERE 
WITH APPLICATION TO EARTH SATELLITES: 
THE SPINNING PLATE 


by 

Fred G. Cunningham 
Goddard Space Flight Center 


INTRODUCTION 

In a previous paper* (Part I) the power input to a stationary flat plate was calculated. The pres- 
ent paper will apply these results to the problem of the spinning plate in order to determine the in- 
stantaneous earth- radiated power incident upon the solar cell paddles associated with the various 
satellites. This analysis does not consider any shielding of the paddles by other paddles or by the 
body of the satellite itself. 

In the previous paper we defined three cases for which the limits of integration were easily de- 
terminable. These were dependent upon the value of the angle between the normal to the plate and 
the earth- satellite line. For a spinning plate, however, this simplification is not possible since a 
spinning plate whose axis of spin is not coincident with the normal to the plate defines a varying an- 
gle between the normal and the earth-satellite line. 

In general, a free rigid body whose axis of spin does not coincide with a principal axis exhibits, 
to an external observer, a precession of the spin axis about the angular momentum vector and a pre- 
cession of the axis of symmetry (in this case the normal) about the spin axis, so that both the axis of 
symmetry and the axis of spin appear to precess about the angular momentum vector. However, 
since we are considering plates attached to larger satellite bodies which usually spin about their 
axis of symmetry (also a principal axis), the spin axis coincides with the angular momentum vector 
and the net result is simply a precession of the normal about the spin axis. The main assumption 
here is that the satellite's spin axis deviates no more than slightly from its axis of symmetry and 
hence, possesses a negligible amount of precessional motion. 


’Cunningham, F. G., "Power Input to a Small Flat Plate from a Diffusely Radiating Sphere with Application to Earth Satellites”, NASA 
Technical Note D-710 (Corrected copy), August 1961. 
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From the aforementioned paper, we can write the general expression for the instantaneous power 
input to a stationary flat plate as 

P = 


2AeaA f 


( cos A cos + s in A sin 0 cos 0) sin <p d0 d0 


(i) 


where: 

A = the generalized emittance of the surface (or or 4 for thermal radiation); 
e = the emissivity of the surface of the radiating sphere; 
a = the absorptivity of the plate; 

A = the elemental vector area of the plate; 

< P = the angle defining the position from the earth-satellite line of the radiating element ds on 
the disk replacing the sphere; 

0 = the azimuthal designation of the element ds; and 

A = the angle between the normal to the plate and the earth- satellite line. 

The upper limit of the 0 integration is given by <p m - sin'Ml/H) where H is the earth-satellite vector, 
and H is given in units of mean earth radii. The upper limit of the 0 integration is determined by the 
conditions governing the three previously mentioned cases: 

(1) 0 - A - (tt/2 ~ 4 > m ) 

(2) (tt/ 2 -0 m ) < A - 77/2 

(3) tt/2 < A < (tt/2+<0 - 

In the first case the upper limit of 0is 77 for all 0, and in the second case the upper limit is v for 
0 - 0 - (tt/2 - a). For the remaining portion of the second case and for all of the third case the upper 
limit of 0 is given by 

“ cos A cos <p /r ,\ 

cos 0 m - sinAsin 0 * 

EXTENSION OF THE METHOD TO A SPINNING PLATE 

The flat plate whose spin axis u passes through the center of mass 
but is not coincident with the normal A is shown in Figure 1 , 
where 



Figure 1 —Geometry showing 
the orientation of the spin 
axis and normal to the plate 


2 


H = the earth-satellite vector; 

a = the angle between A and 

b = the angle between H and «. 

We now define a new angle £ as the azimuthal angle of spin 
of the plate about w ; i.e., w = d£/dt . Figure 1 shows the 
zero value of the spin angle £ for the values of a and b 
shown. In general the zero value of £ is defined when A lies 
in the H,« plane and X is a minimum. Because of symmetry 
we need not consider £ outside the range 0 - £ - tt. In ad- 
dition, w is assumed to be constant and of sufficient magni- 
tude so that we need not consider a time weighted average 
over the angle £ to determine the average power input over 
the spin period, which is our ultimate goal. 


w 



Figure 2— Unit sphere from which the 
angle A. can be calculated 


We shall consider the vectors «, H and A as unit vectors. Figure 2 shows the general orientation 
of the spinning plate, where the ends of the vectors w, H and A lie on the surface of a unit sphere. 

With the above definitions, we can employ spherical trigonometry to determine K Clearly, 


cos A. ~ cos a cos b + sin a sin b cos £ • 


(3) 


By inserting Equation 3 into Equation 1, we have 


p 




cos a cos b 


+ sin a s in b cos £ ) s in <fi cos <$> dp d p 


2\eaA 


JJi- 


i cos a cos 


b + sina s in b cos £) 2 ] 1/2 cos sin 2 <p dp dp 


(4) 


The upper limit of the P integration as given by Equation 2 becomes 


cos p 

•m 


[i 


-(cos a cos b + sin a sin b cos £) cos <p 
( cos a cos b + sinasinb cos £ ) 2 ] 1/2 sin <p 


(5) 


and 


s in p 

f m 


sin 2 4> ~ (cos a cos b + sin a sin b c 
[l - (cos a cos b + sina sinb cos £ ) 2 (] 


wcny 

sin 2 < i>) 


( 6 ) 


However, as pointed out in the earlier paper, there are values of 4> for which the argument of Equa 
tion 5 is greater than unity although X lies within the proper bounds; for this situation 4> m = The 
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values of <p in question are listed on page 2, A complete discussion of the values of 4> and A for which 
Equations 5 and 6 are applicable will be presented subsequently. 

Performing the \p integration of Equation 4 we have: 


2AeaA t 

~ { cos acosb + sinasinb cos £ ) \p sin <p cos <fi d<p 


2Aea-A 


r < P m 

[l — (cos acosb + sinasinb cos s> 2 ] 1/2 sin^ sin 2 <£ d<p 

* n 


(7) 


In the regions where Equations 5 and 6 apply we have, 


C*m 

I (cos a cos b + sin a 

f 

” ( cos a cos b + sin a s in b cos £ ) cos \ 

**0 

s in b cos £) cos ’S ^ 

“ ( cos a cos b + sin a sin b cos £ ) 2 ] 1/2 sin 


s in cos <p d<fi 


+ 


2AeaA 

77 



sin 2 cf> 


(cos acosb + sinasinb cos 


s ) 2 ] 1/2 


sin <fi d<fi 


(8) 


The values of A which define the upper bounds of the regions of interest are \ = 77/2 , 77/2, 

tt/ 2 +<f> m , for which the corresponding £ (for a given a and b) can be determined. These £ are defined 
as follows: 


sin$ m - cos a cos b 

cos C, = sTrTa sin b ■ for X = "/2 - 4> a ; 

v . - cos a cos b 

cos £ ” sin a sin b » ^ ; 

- (s in c/> m + cos a cos b) 

cos ^2 ~ sin a sin b • ^ 0r ^ " 77 + * 

Equation 8 gives the instantaneous power on a spinning flat plate as a function of the altitude 
(from the upper limit of the 4> integration), the angle between the instantaneous spin axis and the nor- 
mal to the plate a, the angle between the instantaneous spin axis and the earth- satellite line b, and 
the azimuthal angle of spin £ about the spin axis w. In the equation the arc cosine factor in the first 
term and the square root factor in the second term arise from the \p integration; the first corre- 
sponds to v/> ro and the second to sin \p m . The arguments of these two factors contain functions of <£, 
and functions of a, b, and £. These latter factors arise from the value of cos k given by Equation 3, 
In the previous paper the value of k was assumed to be constant and calculations were made in the 
appropriate range for various values of A. In that case it was shown that for certain values of A the 
range of integration over <t> had to be divided into two subdivisions, for the first x P m - -rr and the abso- 
lute value of the argument of the arc cosine factor is greater than unity, and for the second, is 


(9a) 

(9b) 

(9c) 
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less than -n and the correct value is given by the arc cosine factor. In fact, since k was held con- 
stant, it was possible to define the limiting value of the 4> integration which depended upon k, and 
hence, to break the integral into its parts — only one of which contained the arc cosine factor as an 
upper limit. However, in the present discussion we are considering the plate as it spins about an 
arbitrarily oriented axis. Consequently, it follows from Equation 3 that Ms no longer constant in 
time but varies with £ for the constant angles a and b. Then, in order to define the limiting values of 
the 4> integration we must make use of the dependence of k upon £ where the corresponding limiting 
values of £ are given by Equation 9. 

As mentioned previously, a plate spinning about an axis not coincident with its normal (which is 
a principal axis) exhibits a wobbly motion as the axis of symmetry (the normal) precesses about the 
axis of spin. Before, we were able to hold k within the range for which the power is incident upon 
the side of the plate in question. The upper limit of this range is given by Equation 9c. When the 
combination of the angles a, b, and £ is such that k, as given by Equation 3, is greater than 77/2 + 4 > m 
then the side of the plate under consideration is not visible from the earth, and the incident power 
is zero. 

Clearly we can define three distinct combinations of the angles a and b: 


(1) 

k $ (77/2 for 0 - £ - 77 ; 

(10a) 

( 2 ) 

0 - x < (77/2 +4> m ) for 0 - £ - £ 2 ; 

(10b) 

( 3 ) 

k > (77/2 +<£ ro ) for 0 - £ - 77 . 

(10c) 


Equation 10c defines the situation with no radiation incident on the side of the plate in question; and 
is of no further interest here. The condition defined in Equation 10a is the least complicated ex- 
pression, integrable over £, for the power because the upper limit of the ^integration is v for all 
values of 4> in the range 0 < <i> m . The general case is represented by Equation 10b where we in- 

tend to imply that, over the range of £ from 0 to n, the value of k might be less than n/ 2-<t> m for cer- 
tain values of £, greater than V 2 but less than tj/2 +4> m for other values of £, and greater than 
7T / J + f° r other values of £, or any other variation on the theme. 

Equation 8 is the expression which must be integrated over <£. It is identical to the equation in 
Part I except that the corresponding expressions involving a, b, and £ for cos k and sin \ are included. 
Before attempting to calculate a given problem we must first examine the range of values of k (given 
by Equation 3) to determine the ranges of £ which are of interest. After the 4> integration it then be- 
comes necessary to average over £, being careful to keep the appropriate range of £ coupled to the 
proper expression. 

For the case where Equation 10a applies, Equation 8 integrates to 


P 


= 2AeaA 


r 


sin"" 1 < 1/H) 

(cos a cos b 


+ sinasinbcos £ ) sin $ cos ^ d<f> 


(ii) 
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which becomes 


p _ ^ €g A , ( cos a cos b + sin a sin b cos £) . (12) 

H 2 

Equation 12 is identical to the expression given as P* in Part I. We indicate the average over £ 
by <P> 4 which is 


<P> 


c 


AeaA 

H 2 


cos a cos b 



(13) 


AeaA 

H 2 


cos a cos 


b 


(14) 


The most general integral expression for the case of Equation 10b is: 


P 


t>» - 1 <i/w) 

t cos a cos b + s in a s in b cos £) sin <p cos <f> &t> 


+ 2AeoA r (cos a cos b ♦ sin a sin b cos 4) sin <£ cos <t> <$£ 



^ pin 1 (l/R) r 

- I ( cos a cos b + sinasinb cos 4 ) cos' 1 < 
Jjt/2 - [coi'* (r« * to. b + tin i tin b cot {>] t 


( cos a cos b + sinasinb cos 4 ) c os <p 
co 5 a cos b ♦ sin a s in b cos C! 3 1' 3 sin 


in 1 (i/K) 

— ^ — I [sin 3 * ( cos a cos b + sinasinb cos 0*] ,/3 sin ^ d<£ 

J n/J - fctt l (toi i n>i b + lir I lid beat f)j 


4 A /•firt ' 1 £l/H) 

2AeaA j , . , . , ,, 

+ — - — j {cos a cos b + sin a s in b cos 4 ' cos 

Jett 3 (cot t tot b * lilt ■ *in b cot C) ” B /J 


‘{l^l 


j»* in ” 1 (I/H> 

i - 2r ^ a ^ I [sin 3 4 s ~ (cos 3 cos b + sl n a sinbeos 4) 2 J V * sin<£ < 

Jcdi ' 1 (cot i e«i b i tin « tinb e»t f) " w/J 


— > sin 4> cos 4> < 


-{cos a cos b + sin a sin b cos 4 j cos# 
cos a cos b + Sin a sinb cos 4) a ] l/ * Hn <j5 


f sin 4> cos 


■4 ' \ 


H ^ . {- cos a cos b\ i 

~J < £ ” cos ^ sin a s in b j f 


_ J ~ cos a cos b \ 
* \ sin a sin b j 


< 4 - cos' 


- cos a cos b “ n 


(15) 


Equation 15 integrates out exactly as Equations 21 and 35 of Part I except that we have replaced the 
cosine k and sine k factors by their equivalents in terms of a, b , and £ . 
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We now have: 


p = 


AeoA 


( cos a cos b + sina sinbcos U |o - £ - £jj 


AeaA 


I - 


_ (" 2 - »)' 


[l - { cos a cos b + sinasinb cos £ ) 


ap7r| 


If . j -(I! 2 “ 1 ) 1 2 (cos a cos b + sin a s in b cos c| 

+ 777 (cos a cos b + sinasinb cos £ ) cos '\~r ~ — — r ; . ~ t . /3 — 

H* \ I [1 - (cos a cos b + sin a s in b cos £ ) 


, . r TJ/2 

(H 2 — l) 1/2 1 - H 2 (cos a cos b + sin a sin b cos £} 2 J 




AeaA 

77 

7 T 

~2 






( H 2 - l ) 1/: * 


HI - (cos a cos b + sin a sinbcos 


_J , 

u 2 T i 77 } + f (HJ _1 ) ,/J [l ~ H 2 ( cos a cos b + sinasinb cos £ } 2 ] l/ 


+ (cos a cos b + sin a sin b cos £) cos 


~(h 2 ~ l) 1/2 (cos a cos b + sinasinb 
[l - (cos a cos b + sin a sin b cos £ ) 


i cos £ )\ ] 

1 J T7T ~J j 


{i* < « 5 O 


(16) 


where the expression in parentheses following certain terms in Equations 15 and 16 indicates the 
range of i over which the preceding terms are applicable. Because of the complexity of some of the 
terms in Equation 16 the average over £ cannot be easily calculated in closed form so that this por- 
tion of the problem will be left to the digital computer. 

Equation 16 applies in its entirety only to the situation for which \does, indeed, have values 
within all three regions defined by o - k - (tt/ 2 ~4> m ), (tt/ 2 -^) < \ - v/2 , and tt /2 < \ < (n/2 +^ m ) each 
of which corresponds to a particular range of values of C as £ goes from 0 to tt. If \ does not take on 
values in all the regions for the entire range of £ certain modifications of Equation 16 must be made 
which are pointed out subsequently. To determine which parts of Equation 16 to use and the corre- 
sponding limits for any case, we must examine the range of \ for that case by checking its extremes. 
Since A precesses about « at a constant angle a, the extremes can be determined by considering the 
cases when A lies in the H, «plane; i.e., when £ = 0 , 77 . The only additional requirement is that k - n. 
In defining values of a and b it suffices that if we let b lie in the range 0 - b - 77 , we need only con- 
sider a in the range 0 - a 5 77/2. For example, the situation represented by a - 100°, b = 10° is 
identically represented by the situation b = 170°, a = 80°. The direction of «, and consequently the 
value of b, is always determined by the right hand rule. It also follows that interchanging the values 
of a and b does not alter the problem. In addition, if we wish to determine the input to the other side 
of the plate (called side /3) after having determined it for one side (called side a), by the methods 
outlined above, the extremes of the range of k for side p can be easily calculated by subtracting the 
limits for side a from 180 degrees. Then, once the extremes of k fi are known, it is quite easy to 
construct the corresponding values of a and b which, used in conjunction with Equation 3, will allow 
for the determination of all the parameters of the new problem. For example, if a = 80°, b = 170°, 
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for side a, we have that 90° - A - 110°; whereas 70° 5 X - 90° for side p. The latter set corresponds 
to b = 10°, a = 80° or a = 10°, b = 80°. 

We now outline the various possibilities of the general case represented by Equation 10b. But 
first we shall make the following definitions for simplicity: 


A 


AeaA 

H 2 "" 


f cos a cos b 


+ 


sin a sin b cos 


c] . 


and 


AeaA 

77 _ , 

j 

(h 2 - 

i) i/2 __i 

77 

■y ” sin 

\h[i - 

(cos a cos b + 

sinasinb cos £ ) 2 J 1/2 j 


H 2 


(cos a cos b + sin a sin b cos £) cos 


_ ^ J" - |^H 2 ~ l) 1/2 ( cos acosb + sin a sinb cos J] 
"p - {cos acosb + sin a sinb cos £ ) 2 J 1/2 


-(h 2 -l) l/a [l 


H 2 (cos acosb + sin a sinbcos 



Equation 16 can now be written in symbolic form: 

P = a(o - $ - £,) + b(c, < £ - c) + b(c' < £ - Q ■ 

The first step in any problem is to determine the range of A. We shall define three cases: 
A - 77/2 for all £ ; 

A 7 77/2 for suitable values of £ in the range 0 - £ - 77; 

A > 77/2 for all £. 


(17a) 


(17b) 


(18) 

(19a) 

(19b) 

(19c) 


In 19a the only value of £ of concern is that for which A = tt/2 -<f> m or £ = £, as given by Equation 9a. 
If for a given value of H the argument of the arc cosine in 9a is greater than +1 then A > ( 77/2 ~4> m ) 
for all £ and Equation 16 reduces to 


P = B ( 0 - £ - 77) • ( 20 ) 

If, for a given value of H, the argument of the arc cosine in Equation 9a is greater than or equal to -1 
but less than or equal to +1, so that Equation 9a defines the value of £ at which A = (77/2 then 

Equation 16 becomes 


P = A (0 - £ - £j + b(£j < £ - 77) • (21) 

However, if for a given value of H the argument of the arc cosine term in 9a is less than 1, then 
A < (77/2 for a H £ Equation 16 becomes 
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p 


A(0 - £ - 7 t ) , 


( 22 ) 


which, when averaged over C, yields Equation 14 once again. 

From 19b it immediately follows that X can never be less than (V 2 or greater than (n/2 +4> m ) 
for all £. Equations 9a and 9c define the limits of £ in question. In the preceding statement, the fact 
that the value of the argument of the arc cosine in 9a can never be less than -1 and that the value of 
the argument of the arc cosine in 9c can never be greater than +1 can easily be shown for any prob- 
lem satisfying the conditions of 19b. If the argument of the arc cosine in 9a is greater than or equal 
to —1 but less than or equal to +1 and the argument of the arc cosine in 9c is less than —1, Equa- 
tion 16 reduces to 

P = A(0 < t * O + B(Cj < i < J . (23) 

If both the arguments of the arc cosines in 9a and 9c lie between +1 and -1 we then have Equation 16. 
If the argument of the arc cosine in 9a is greater than +1 while the argument of the arc cosine of 9c 
lies between +1 and — 1, Equation 16 becomes 

P = B (o - £ - £ 2 ) ( 24 ) 

Fi nall y, if the argument of the arc cosine in 9a is greater than +1 while the argument of the arc co- 
sine in 9c is less than -1, Equation 16 reduces to 

p = B( 0 - £ - n) • (25) 

In 19c the only value of £ of concern is that for which X = rr/2 +4> n or £ = £ 2 as given by Equation 9c. 

If the argument of the arc cosine in 9c is less than -1, x< (V2 + <t> m ) for all £ and Equation 16 
becomes 


P = B(0 - £ - 7 r) ■ (26) 

If the argument of the arc cosine of 9c lies between +1 and -1, the value of £ 2 is defined and Equa- 
tion 16 is 


p = b(o - £ - £ 2 ) ■ (27) 

If the argument of the arc cosine of 9c is greater than +1 then X > (n/2 + <£ m ) for all £, and the corre- 
sponding power input is zero. 

DISCUSSION 

Equation 16 is general and applies to any sort of diffusely radiating sphere. However, in the fol- 
lowing calculations we are considering the diffusely radiating sphere to be the earth which radiates 
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as a blackbody at a temperature approximately equal to 25(TK. Thus 


where e equals unity because of the equivalent blackbody approximation and a is taken equal to unity 
in order to give the energy incident upon the plate. In addition, Equation 16 can be modified to give 
directly the incident power to the plate as a function of the orbital position of the satellite by giving 
H as a function of the orbital parameters and the azimuthal position of the satellite in orbit. This 
can be done by utilizing the polar equation of an ellipse, 


_ a (l - e 2 ) 

1 + e cos 77 ’ 


(29) 


where n is the angle giving the azimuthal position of the satellite from perigee, a is the semi- major 
axis and e is the eccentricity. However, in the subsequent evaluations of Equation 16 this additional 
information will not be included because the presence of the arc cosine factor in the equation renders 
it almost impossible to write, in terms of the orbital parameters, the equation in closed form (after 
averaging over 5 ) so that not a great deal is to be gained unless we wish the computer to print out 
the value for each orbital position. 

So far we have only considered the case where the spin axis passed through the center of the 
plate. However, it is easy to see that the dependence of k on £ remains identical to that given by 
Equation 3 for the case where a and b are defined as usual but the axis of spin is displaced from the 
center of the plate. As long as this lateral displacement (usually half the width of a satellite) is 
small in comparison to the distance of the plate from the surface of the sphere (which always ob- 
tains for earth satellites) the previous analysis remains applicable. 


RESULTS 

In the graphs that follow, the average incident power to the Hat plate per spin period is plotted as 
a function of altitude above the surface of the earth for various values of the angles a and b. The situ- 
ation occurring whenb equals 0 or n results in \ = a for all 4 so that these values are not considered, 
since the incident power for \= constant is readily available in Part I. This situation also applies 
for a = 0 for which \ = b = constant. The first 16 graphs apply for H in the range 200 5 h 5 3200 km 
(Appendix A), while the remaining graphs apply to the entire range 200 $ h -32000 km (Appendix B). 

ACKNOWLEDGMENT 

The author wishes to thank Mr. E. Monasterski, Goddard Space Flight Center, for the IBM 7090 
computations. 


10 



Appendix A 

The Average Incident Power to a Flat Plate per Spin Period for 
Various Angles a and b as a Function of Altitude for the 
Range 200 4 H l 3200 km 
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Figure A16-b = 160°, 170' 






Appendix B 

The Average Incident Power to a Flat Plate per Spin Period for 
Various Angles a and b as a Function of Altitude for the 
Range 200 4 H 4 32000 km 
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Figure B1 — b = 10 



